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Abstract: A search for high-mass resonances decaying into final states using 

proton-proton collisions at ^/s = 8 TeV produced by the Large Hadron Collider is pre¬ 
sented. The data were recorded with the ATLAS detector and correspond to an integrated 
luminosity of 19.5-20.3 fb“^. No statistically significant excess above the Standard Model 
expectation is observed; 95% credibility upper limits are set on the cross section times 
branching fraction of Z' resonances decaying into pairs as a function of the reso¬ 

nance mass. As a result, Z' bosons of the Sequential Standard Model with masses less 
than 2.02 TeV are excluded at 95% credibility. The impact of the fermionic couplings on 
the Z' acceptance is investigated and limits are also placed on a Z' model that exhibits 
enhanced couplings to third-generation fermions. 
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1 Introduction 

Searches for new heavy resonances decaying to tau lepton pairs are both theoretically and 
experimentally well motivated [1-6]. Heavy Z' bosons often arise in grand unified theories 
and while they are typically considered to obey lepton universality, this is not necessarily 
a requirement. In particular, some models offering an explanation for the high mass of the 
top quark predict that such bosons preferentially couple to third-generation fermions [7, 8]. 
Models containing non-universal Z' bosons can explain the anomalous dimuon production 
observed at the DO experiment [9, 10] and the excess in semileptonic R-meson decays into 
tau leptons observed at the Belle and BaBar experiments [11-13]. Searches in the ditau 
channel are also sensitive to sgoldstino-like scalars in supersymmetric models [14, 15], 
hidden sector Z' models [16] and to the anomalous tau lepton dipole moments and higher- 
order tau-gluon couplings [17]. 

In this article, a search for high-mass resonances decaying into final states using 

proton-proton (pp) collisions at a center-of-mass energy of -y/s = 8 TeV produced by the 
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Large Hadron Collider (LHC) [18] is presented. The data were recorded with the ATLAS 
detector [19] and correspond to an integrated luminosity of 19.5-20.3 fb“^. Tan leptons 
can decay into a charged lepton and two neutrinos (riep = Te or r^), or hadronically (rhad), 
predominantly into one or three charged pions, a neutrino and often additional neutral 
pions. The Thad'^hadi '^M'^had and TeThad channels are analysed, accounting for 42%, 23% 
and 23% of the total branching fraction, respectively. A counting experiment is 

performed in each channel from events that pass a high-transverse-mass requirement. Due 
to the different dominant background contributions and signal sensitivities, each channel 
is analysed separately and a statistical combination is used to maximise the sensitivity. 

The Sequential Standard Model (SSM), which contains a boson with couplings 

identical to the Standard Model Z boson, is chosen as the benchmark model to optimise 
the analysis and to quantify the experimental sensitivity. Limits on the -^gg^ cross section 
times the branching fraction in tau pairs, cr{pp — >■ -^gg^ + X) ■ H(Zggj^ —>■ T+r“) = crBssM, 
are provided as a function of the resonance mass, mz'- The impact on the signal acceptance 
times efficiency from changing the .^gg^ couplings is assessed, which allows the limits on 
ZggM to be reinterpreted for a broad range of models. Limits are also placed on the non- 
universal G(221) model [8, 20, 21], which contains a Z'-^^ boson that can exhibit enhanced 
couplings to tau leptons. 

Direct searches for high-mass ditau resonances have been performed by the ATLAS 
and CMS collaborations using 5 fb“^ of integrated luminosity at ^/s = 7 TeV [22, 23]. 
The searches exclude .^gg^ with masses below 1.4 TeV at 95% CL.^ For comparison, the 
most stringent limits on ^ggjy; in the dielectron and dimuon decay channels combined 
are 2.90 TeV at 95% CL from both ATLAS [24] and CMS [25]. While the limits on 
a{pp — )■ .^ggjvi + X) ■ B{Z' — )■ eXe~jpXpT) are in general stronger than those on crHggM, 
they may be evaded by models with weak couplings to electrons and muons. Indirect lim¬ 
its on Z' bosons with non-universal flavour couplings have been set using measurements 
from LEP and LEP H [26] and translate to a lower bound on the Z' mass of 1.09 TeV at 
95% CL. Indirect limits have also been placed on the non-universal G(221) model [8, 27-29]. 
The strongest exclude Z'-^-^j with a mass lower than 1.8 TeV at 95% CL. 

This article is structured as follows. Section 2 provides an overview of the ATLAS 
detector. The event samples used in the analysis, recorded by the ATLAS detector or 
simulated using the ATLAS simulation framework, are described in section 3. The recon¬ 
struction of physics objects within the event samples is described in section 4. A description 
of the selection criteria used to define Z' signal regions is given in section 5. Section 6 de¬ 
scribes the estimation of background contributions, followed by a description of systematic 
uncertainties in section 7. In section 8, the impact of altering the Z' couplings on the signal 
acceptance is described and the non-universal G(221) model is introduced. A presentation 
of the results is given in section 9, followed by concluding remarks in section 10. 


^CL is used interchangeably throughout this article to refer to both confidence level (frequentist) and 
credibility limit (Bayesian). 
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2 ATLAS detector 


The ATLAS detector at the LHC covers nearly the entire solid angle around the colli¬ 
sion point. It consists of an inner tracking detector surrounded by a thin superconducting 
solenoid, electromagnetic (EM) and hadronic calorimeters, and a muon spectrometer in¬ 
corporating large superconducting toroid magnets. 

The inner-detector system is immersed in a 2 T axial magnetic field and provides 
charged-particle tracking in the range \'q\ < 2.5.^ A high-granularity silicon pixel detector 
covers the vertex region and typically provides three measurements per track. It is followed 
by a silicon microstrip tracker, which usually provides four pairs of measurements per track. 
These silicon detectors are complemented by a transition radiation tracker (TRT), which 
enables radially extended track reconstruction up to jr/l = 2.0. The TRT also provides 
electron/pion discrimination based on the fraction of hits (typically 30 in total) above a 
higher energy-deposit threshold corresponding to transition radiation. 

The calorimeter system covers the pseudorapidity range \ri\ < 4.9. Within the region 
\r]\ < 3.2, EM calorimetry is provided by high-granularity barrel and endcap liquid-argon 
(LAr) EM calorimeters with lead absorbers, with an additional thin LAr presampler cov¬ 
ering \r]\ < 1.8 to correct for upstream energy loss. Hadronic calorimetry is provided by a 
steel/scintillator-tile calorimeter, segmented into three barrel structures within \r]\ < 1.7, 
and two copper/LAr hadronic endcap calorimeters. Coverage in the forward region is 
achieved by copper/LAr and tungsten/LAr calorimeter modules optimised for EM and 
hadronic measurements, respectively. 

The muon spectrometer comprises separate trigger and high-precision tracking cham¬ 
bers measuring the deflection of muons in a magnetic field generated by superconducting 
air-core toroids. The precision chamber system covers the region \r]\ < 2.7 with three layers 
of monitored drift tubes, complemented by cathode strip chambers in the forward region, 
where the background is highest. The muon trigger system covers the range \r]\ < 2.4 with 
resistive plate chambers in the barrel, and thin gap chambers in the endcap regions. 

A three-level trigger system is used to select interesting events [30]. The Level-1 trigger 
is implemented in hardware and uses a subset of detector information to reduce the event 
rate to a design value of at most 75 kHz. This is followed by two software-based trigger 
levels which together reduce the event rate to a maximum of 1 kHz. 

3 Event samples 

The data used in this search were recorded with the ATLAS detector in pp collisions at a 
centre-of-mass energy of ^/s = 8 TeV during the 2012 run of the LHC. Only data taken 
with pp collisions in stable beam conditions and with all ATLAS subsystems operational 
are used, resulting in an integrated luminosity of 20.3 fb“^. Eor the analysis of the ThadRiad 

^ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the 2 -axis along the beam pipe. The a:-axis points from the IP to the centre of 
the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, 0) are used in the transverse plane, 
(j> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle 
6 as r] = — lntan(6/2). The geometrical distance between objects is defined as AR = y^(A(^pT-~(A? 7 E. 
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channel, a small fraction of data from the initial running period are discarded as the trigger 
conditions are not accounted for by the simulation, resulting in an integrated luminosity of 
19.5 fb“^. The ThadThad channel uses events passing a single-tau trigger with a transverse 
momentum (pt) threshold of 125 GeV, designed to select hadronic tau decays. The r^Thad 
channel uses events passing a single-muon trigger, either with a threshold of 24 GeV 
including an isolation requirement or with a threshold of 36 GeV without an isolation 
requirement. The TeThad channel uses events passing a single-electron trigger, either with a 
Pt threshold of 24 GeV including an isolation requirement, or with a threshold of 60 GeV 
without an isolation requirement. Events that pass the trigger are selected if they contain 
a vertex with at least four associated tracks, each with px > 0.5 GeV. Events may have 
several vertices satisfying this requirement due to multiple pp interactions occurring in the 
same or neighbouring bunch crossings, referred to as pile-up. The event vertex is chosen 
as the one with the largest sum of the squared track transverse momenta. 

Monte Garlo (MG) simulation is used to estimate signal efficiencies and some back¬ 
ground contributions. Simulated samples of events from the following background processes 
are used: Zj^* —)■ rr and Z/ 7 *(—)■ .^.^)-|-jets {I = e,//) enriched in high-mass events, and 
VE-|-jets, tt, single-top-quark {Wt, s-channel and t-channel) and diboson {WW, WZ, and 
ZZ) production. Each sample is produced with one of the following event generators: 
Pythias. 165 [31], Sherpa 1.4.1 [32], MC@NLO4.01 [33-35], AcerMC3.8 [36], Her- 
wiG 6.520 [37] or PowHeg-Box 1.0 [38-41]. The most consistent set of available samples 
was chosen. The Z/y* —>■ rr process is generated at leading order so that the sample can 
also be reweighted to describe the Z' signal. The combination of ti and single-top-quark 
production are referred to as top. In some cases the generators are interfaced to the follow¬ 
ing external software for parton showering, hadronisation and multiple parton interactions: 
Pythias, Pythia 6.421 [42] or Herwig (which is itself interfaced to Jimmy 4.31 [43] for 
multiple parton interactions). The tau lepton decay is performed by either Pythia 8, 
Sherpa or Tauola [44]. For Pythia 8, the sophisticated tau decay option is used, which 
provides fully modelled hadronic currents with spin correlations for tau-lepton decays [45]. 
In all samples other than those generated with Sherpa, final-state photon radiation is per¬ 
formed by Photos [46]. The CTEQ6L1 [47] and CTIO [48] parton distribution functions 
(PDFs) and the AU2, AUET2, AUET2B [49] and CTIO [32] MG tunes are used. A 
summary is given in table 1. 

The contributions from simulated processes are normalised using theoretical cross sec¬ 
tions. The Z/y* cross section is calculated up to next-to-next-to-leading order (NNLO) in 
QCD including next-to-leading order (NLO) electroweak corrections using FEWZ 3.1 [50] 
configured with the MSTW2008NNLO PDF set [51]. This cross section is used to derive 
mass-dependent A-factors that are used to weight the simulated Z/y* samples. Cross 
sections for the other background processes are calculated without the use of differential 
A-factors to at least NLO in QCD, as specified in table 1. 

The contributions of the various Z' signal models are estimated by reweighting the 
Z/y* —>■ rr sample using TauSpinner [60-62], which correctly accounts for spin effects 
in the tau decays. The algorithm relies on a leading order approximation in which spin 
amplitudes are used to calculate the spin density matrices for hard 2 —)• 2 Born level 
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Process 

Generator 

PS-bMPI 

Tau decay 

PDF set 

MC tune 

Cross section 

Z/y* —>■ TT 

Pythia 8 

Pythia 8 

Pythia 8 

CTEQ6L1 

AU2 

NNLO [50] 

IF-l-jets 

Sherpa 

Sherpa 

Sherpa 

CTIO 

CTIO 

NNLO [52, 53] 

ti 

Single top 

MC@NLO 

Herwig 

Tauola 

CTIO 

AUET2 

~NNLO [54-56] 

(Wt) 

MC@NLO 

Herwig 

Tauola 

CTIO 

AUET2 

~NNLO [57] 

(s-channel) 

MC@NLO 

Herwig 

Tauola 

CTIO 

AUET2 

NNLL [58] 

(t-channel) 

AcerMC 

Pythia 6 

Tauola 

CTEQ6L1 

AUET2B 

~NNLO [57] 

Diboson 

Herwig 

Herwig 

Tauola 

CTEQ6L1 

AUET2 

NLO [59] 

Zh* U 

PowHeg-Box 

Pythia 8 

Pythia 8 

CTIO 

AU2 

NNLO [50] 


Table 1. Details regarding the MC simulated samples. The following information is provided 
for each sample: the generator of the hard interaction, the parton shower and hadronisation (PS), 
multiple parton interactions (MPI) and the tau decay; the PDF set; the MC tune and the order 
in QCD of the cross section calculation. All cross sections are calculated at either NLO, NNLO, 
approximate NNLO (^NNLO) or next-to-next-to-leading logarithm (NNLL). 


processes. The impact of interference between Z' and Zj^* is typically small (as discussed 
in section 8.1), so it is not included. For each signal model, several mass hypotheses are 
considered, ranging from 500 to 2500 GeV in steps of 125 GeV. 

All generated events are propagated through a detailed Geant4 simulation [63] of the 
ATLAS detector and subdetector-specihc digitisation algorithms [64] and are reconstructed 
with the same algorithms as the data. Pile-up is simulated by overlaying minimum-bias 
interactions generated with Pythia 8 (with an MC tune specific to the LHC [65]) on the 
generated signal and background events. The resulting events are reweighted so that the 
distribution of the number of minimum-bias interactions per bunch crossing agrees with 
data. Due to the high momenta of the tau decay products, however, pile-up has little 
impact on the analysis. The effective luminosity of most simulated samples is at least as 
large as the integrated luminosity of the data; the statistical uncertainty from the limited 
sample size is accounted for in the statistical analysis. 

4 Physics objects 

In this section the reconstruction of electrons, muons, hadronic tau decays and the missing 
transverse momentum is described. Preliminary selections are applied to all electrons, 
muons and tau candidates. Further selection is applied to some of the objects as part of 
the event selection described in section 5. Gorrections are applied to the kinematics and 
efficiencies of reconstructed electrons, muons and hadronic tau decays in simulated samples 
so that they match the performance measured from the data. 

The reconstruction, energy calibration and identification of hadronic tau decays in 
ATLAS is described in detail in ref. [66]. Candidates for hadronic tau decays are built 
from jets reconstructed using the anii-kt algorithm [67, 68] with a radius parameter value 
of 0.4. The jets are calibrated to the hadronic energy scale with correction factors based 
on simulation and validated using test-beam and collision data [69] . Only the visible tau- 
decay products (all products excluding neutrinos), Thad-visj are considered when calculating 
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kinematic properties. The calculation of the four-momentum uses clusters with AR < 0.2 
from the initial jet-axis and includes a final tan-specific calibration derived from simulated 
samples, which accounts for out-of-cone energy, energy lost in dead material, underlying- 
event and pile-up contributions and the typical composition of hadrons in hadronic tan 
decays. The size of the tan-specific calibration is typically a few percent. The calibrated 
energy scale in data and simulation have been compared and agree within the ~1.5% uncer¬ 
tainty of the measurement. Candidates are required to have either one or three associated 
tracks (prongs) reconstructed in the inner detector. The tan charge is reconstructed from 
the sum of the charges of the associated tracks and is required to be ±1. The charge 
misidentification probability is found to be negligible. Hadronic tan decays are identified 
with a multivariate algorithm that employs boosted decision trees (BDTs) to discrimi¬ 
nate against quark- and gluon-initiated jets using shower shape and tracking information. 
Working points with a tan identification (ID) efficiency for l-prong/3-prong candidates of 
about 55%/40% (medium) for the r^Thad and TeThad channels and 65%/45% (loose) for the 
ThadWad channel are chosen, leading to rates of false identification for quark- and gluon- 
initiated jets of below a percent. The tan ID efficiency is independent of pT and pile-up. 
Corrections of a few percent are applied to the efficiency in simulation. Candidates arising 
from the misidentification of electrons are rejected using a separate BDT. In the r^Thad 
channel, a dedicated selection is applied to suppress candidates arising from the misiden¬ 
tification of muons. Tan candidates are required to have pT > 30 GeV and to be in the 
fiducial volume of the inner detector, \7]\ < 2.47. The transition region between the barrel 
and endcap EM calorimeters, with 1.37 < |p| < 1.52, is excluded. In the riepThad channels, 
candidates that have the highest-pT track in the range \rj\ < 0.05 are rejected. This region 
corresponds to a gap in the TRT, which reduces the power of electron/pion discrimination. 

Muon candidates are reconstructed by combining an inner-detector track with a track 
from the muon spectrometer. The candidates are required to have px > 10 GeV and 
IpI < 2.5. Muon quality criteria are applied to achieve a precise measurement of the muon 
momentum and reduce the misidentification rate [70]. These quality requirements corre¬ 
spond to a muon reconstruction and identification efficiency greater than 95%. 

Electrons are reconstructed by matching clustered energy deposits in the EM calorime¬ 
ter to tracks reconstructed in the inner detector [71]. The tracks are then refitted using the 
Gaussian Sum Filter algorithm [72], which accounts for energy loss through bremsstrahlung. 
The electron candidates are required to have px > 15 GeV and to be within the fiducial 
volume of the inner detector, \r]\ < 2.47 (the EM calorimeter transition region is excluded). 
The candidates are required to satisfy quality criteria based on the expected calorimeter 
shower shape and amount of radiation in the TRT. These quality requirements correspond 
to an electron identification efficiency of approximately 95% [73] . 

Electrons and muons are considered isolated if they are away from large deposits of 
energy in the calorimeter and tracks with large px consistent with originating from the same 
vertex. Lepton isolation is defined using the sum of the transverse energy, Ex, deposited 
in calorimeter cells with AR < 0.2 from the lepton, and the scalar sum of the px of 

tracks with px > 0.5 GeV consistent with the same vertex as the lepton and with AR < 0.3, 
pSp^. Muons are considered isolated if they have Px^/px < 5%. Isolated electrons must 
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have < 5% and < 5 GeV + 0.7% x px and must pass a tighter identification 

requirement corresponding to an efficiency of approximately 70%. 

Geometric overlap of objects with AT? < 0.2 is resolved by selecting only one of the 
overlapping objects in the following order of priority: muons, electrons, tau candidates. The 
order is determined by the ability to identify the objects from their detector signatures. 
The missing transverse momentum, with magnitude is calculated from the vector 

sum of the transverse momenta of all high-px objects reconstructed in the event, as well as 
a term for the remaining activity in the calorimeter [74]. Clusters associated with electrons, 
hadronic tau decays and jets are calibrated separately. The remaining clusters are weighted 
using tracking information to reduce the effect of pile-up on the 77™'®® resolution. A single 
weight is calculated for each event using all tracks that are not matched to high-px objects. 
The tracks are categorised based on whether or not they are matched to the primary vertex. 
The weight is then defined as the ratio of the sum of the px of tracks originating from the 
primary vertex to the sum of the px of all tracks. 

5 Event selection 

Selected events in the Thad^had channel must contain no electrons with px > 15 GeV or 
muons with px > 10 GeV and at least two tau candidates: one with px > 150 GeV that 
is matched to the object that passed the trigger and the other with px > 50 GeV. This 
constitutes the preselection. If multiple tau candidates are selected, the two highest-px 
candidates are chosen. This decision is made before applying the BDT tau ID, to avoid 
kinematic biases in control regions defined by reversing the ID requirement. The tau can¬ 
didates are then required to have charges of opposite sign (OS). Finally, the angle between 
the tau candidates in the transverse plane, A(/)(ri, r 2 ), must be greater than 2.7radians, as 
tau leptons from the decay of heavy neutral resonances are typically produced back-to-back 
in the transverse plane. 

Selected events in the riepThad channels must contain exactly one isolated muon 
with Px > 30 GeV or one isolated electron with px > 30 GeV; no additional electrons 
with Px > 15 GeV or muons with px > 10 GeV; and at least one tau candidate with 
Px > 30 GeV. This constitutes the preselection. If multiple tau candidates are selected, 
the highest-px candidate is chosen. As in the ThadAiad channel, this choice is made be¬ 
fore applying the BDT tau ID. The angle between the lepton and tau candidate in the 
transverse plane, A(/)(7, r), must be greater than 2.7radians, and they must have opposite 
charge. The transverse mass is defined as: 

mTip'^,p^) = ^J2p^p^ (1 - cos A(/>(p^,p^)), 

where p^ and p^ are two reconstructed physics objects with transverse momenta p^^ and p.^, 
respectively, which subtend an angle of A(j){p^,p^) in the transverse plane. The kF-l-jets 
background is suppressed by requiring the transverse mass of the lepton-77™'®® system, 
mx('7,77™'®®), to be less than 50 GeV. 

The search in all channels is performed by counting events in signal regions with total 
transverse mass above thresholds optimised separately for each signal mass hypothesis in 
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each channel to give the best expected exclusion limits. The same thresholds are found to 
be optimal for all channels. The total transverse mass, is defined as 


"iT*(n, T2, ^ m^(Ti,T2) + m|(ri, + m^(r2, , 

where n and T 2 denote the reconstructed visible decay products of the two tau leptons (e, 
fj, or Tha(j_vis). 


6 Background estimation 

The dominant background process in the Thad'^had channel at high mass is Zj’^* —)• rr, which 
is estimated using simulation. The modelling of the pp —)• Zj^* process has been shown to 
be very reliable by using decays to electrons and muons [24, 75]. Additional uncertainties 
related to the modelling of high-px tau decays are also considered, as described in section 7. 
Multijet production makes a large contribution at low mass in the Thad'T'had channel and 
is estimated by weighting events in data where the subleading tau candidate fails tau ID, 
with fake-factors that parameterise the rate for jets to pass tau ID (section 6.1). Due to the 
relatively large size of the sample that fails tau ID, this procedure provides high statistical 
precision, which is particularly crucial in the high-mass tail. The fake-factors are measured 
from data in a separate control region. Diboson, IT-|-jets, ti, Zl^*{-^ £^)-|-jets, and single- 
top-quark production make minor contributions and are estimated using simulation. To 
improve the modelling of these background processes, events in the simulation that contain 
jets misidentified as hadronic tau decays are weighted by fake-rates measured in a IT-|-jets 
control region in data (section 6.3). 

The dominant background contributions in the riepThad channels come from Z/y* —)• rr, 
which is estimated using simulation, and from processes in which a jet is misidentified as 
a hadronic tau decay. The latter is mainly composed of IT-|-jets events and is estimated 
using fake-factors to weight events in data where the tau candidate fails ID, similarly to 
the procedure in the Thad'^had channel (section 6.2). Diboson, tt, Z/ 7 *(—)■ £t’)-|-jets and 
single-top-quark production in which the tau candidate does not originate from a jet make 
minor contributions and are estimated using simulation. In the following subsections, the 
data-driven background estimates are described in more detail. 

6.1 Multijet background in the ThadDiad channel 

To estimate the multijet background in the Thad^had channel, two control regions are used. 
Events in the first control region are required to pass the same selection as in the analysis, 
except for the subleading tau candidate, which is required to fail the BDT tau ID. The 
multijet contribution in the signal region is estimated by weighting these events with a 
tau ID fake-factor. The fake-factor, /tau-iD; is defined as the ratio of the number of tau 
candidates that pass the BDT tau ID, jypass tau-iD ^ number that fail, tau-iD_ 

The fake-factors are calculated from a second control region that is highly pure in multijet 
events, the multijet control region (described below), and they depend on the pT and track 
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multiplicity, Mrackj of the subleading tau candidate: 


/tau-ID {PT ) -^track) 


jS^pass tau-ID 

iVfail*a-ID(pT,iVtrack) 



multijet 


The fake-factors have no significant dependence on r]. The number of multijet events in a 
bin of pt, A^track and any additional variable that is uncorrelated to the BDT tau ID, x, is 
given by: 


A^multijet (PT) A^track) 3:') /tau-ID(PT) A^track) ^ -^data (PTj A^track; • 

The multijet control region is designed to be as similar to the signal region as possible, 
while avoiding contamination from hadronic tau decays. This is achieved by loosening the 
tau ID requirements. Specifically, the selection for this control region is the same as for the 
signal region except with the following alterations. The BDT tau ID is not applied to either 
tau candidate. Instead of using the single-tau trigger, events are selected using single-jet 
triggers with thresholds ranging from 45 to 360 GeV, with all but the highest threshold 
trigger being prescaled. The px of the subleading tau candidate must be at least 40% of 
the pt of the leading tau candidate (px-balance > 0.4) to avoid bias at low px due to the 
disproportionate fraction of events coming from the unprescaled 360 GeV jet trigger. The 
opposite-sign requirement on the charges of the two tau candidates is removed to increase 
the sample size. 

Figures l(left) and 1 (right) show the fake-factors for 1-prong and 3-prong candidates, 
respectively. Use of these fake-factors relies on the assumption that they are insensitive 
to the alteration of the selection between the signal region and multijet control region. 
Systematic uncertainties on the fake-factors are derived by altering the selection on the 
px-balance, the charge product, and the identification of the leading-px tau candidate. 
These variations modify the fractional contribution of quark- and gluon-initiated jets in 
the sample, leading to large variations in the fake-factors at low px where the composition 
is mixed and little variation to the fake-factors at high px where the sample is quark 
dominated. 

6.2 Jet background in the "nepThad channel 

The background contributions originating from quark- and gluon-initiated jets that are 
misidentified as hadronic tau decays in the riepThad channels are modelled using a fake- 
factor method, similar to that used in the Thad'^had channel. In contrast to the Thad'^had 
channel, the background is dominated by IU-|-jets production, with a minor contribution 
from multijet production. To reduce the sensitivity to the differing fake-factors in W -|-jets 
and multijet events (due to a different qnark/glnon fraction), events failing a very loose level 
of BDT tau identification (corresponding to efficiencies of 98% and 90% for 1-prong and 
3-prong hadronic tau decays, respectively) are rejected. This significantly suppresses the 
gluon contribution, which typically consists of wider jets with higher hadron multiplicity 
which are more readily rejected by the tau ID. In the riepThad channels, there is also a non- 
negligible contribution to the first control region (fail-ID control region) from background 
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Figure 1. Tau ID fake-factors for (left) 1-prong and (right) 3-prong tan candidates, measured in 
the multijet control region of the ThadThad channel. The statistical and systematic uncertainties are 
shown, successively added in quadrature. 


processes containing hadronic tau decays, which is subtracted using simulation. The fake- 
factors are measured in a high-purity VF-l-jets control region and they depend on pxj V 
and Altrack of the tau candidate. The IT-l-jets control region uses the same selection as 
the signal region but with the medium BDT tau ID replaced by very loose BDT tau ID 
and the mx requirement replaced by 70 GeV < mx < 200 GeV. A second control region 
enriched in multijet events is defined, which has a higher fraction of gluon-initiated jets 
and represent an extreme variation in the jet composition. This control region uses the 
same selection as the IT-|-jets control region but the lepton is required to fail isolation, 
the mx requirement is replaced by mx < 30 GeV and < 30 GeV is required. A 30% 
systematic uncertainty is derived from the difference in the fake-factors in the multijet and 
IT-|-jets control regions. Figures 2(left) and 2(right) show the fake-factors measured in 
each of the two control regions in the r^Thad channel, integrated across all \r]\ regions, for 
1-prong and 3-prong candidates, respectively. The fake-factors in the TeThad channel are 
similar. 

Finally, in the riepThad channel, two additional steps are taken to ensure is mod¬ 
elled well by the fake-factor estimate. Firstly, the standard reconstruction treats 

the selected tau candidate in the fail-ID control region as a jet rather than a hadronic 
tau decay. Therefore, the FI™®® is recalculated in the fail-ID control region using the tau 
hypothesis for the selected tau candidate. Following this, a slight bias in the shape of the 
Fjmiss digtribution is corrected for by reweighting in bins of the Fl™^®® projected along the 
direction of the tau candidate. An additional 20% uncertainty is applied to the estimate 
of the jet background event yield obtained after the full event selection, derived from the 
difference in the estimate between applying and not applying the Fi™^®® reweighting. 
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Figure 2. Tau ID fake-factors for (left) 1-prong and (right) 3-prong tan candidates, measured in 
the IF-1-jets control region of the r^Thad channel, integrated across all |ry| regions. The statistical and 
systematic uncertainties are shown, successively added in quadrature. The fake-factors measured 
in the alternative multijet control region are overlaid. 


6.3 Jet background other than multijet in the ThadThad channel 

In the ThadThad channel, backgrounds originating from quark- and gluon-initiated jets that 
are misidentified as hadronic tau decays in processes other than multijet production are 
estimated using simulation (predominantly VF-l-jets). Rather than applying the tau ID 
to the simulated jets, they are weighted by fake-rates. This not only ensures the correct 
fake-rate, but enhances the statistical precision of the estimate, as events failing the tau 
ID are not discarded. The fake-rate for the sub-leading tau candidate, Rtau-ro'^> defined 
as the ratio of the number of tau candidates that pass tau ID, , to the total 

number of tau candidates, The fake-rate for the leading tau candidate, is 

defined as the ratio of the number of tau candidates that pass tau ID and the single-tau 
trigger requirement, to The fake-rates are calculated from a 

second control region that is high in IF-|-jets purity (described below), and they depend 
on pt and -iVtrack of the tau candidate: 

jYpass tau-ID -h trigger Mtrack) 


^llu4D(PT,iVtrack) = 




VF-l-jets 


^t“(PT,iVtrack) = 


_ iVP-^*""-™(pT,iVtrack) 


All simulated events are assigned a weight: 


IT-l-jets 


WMC= n (1 - [1 - ^Lu-Id(P^> ^^rack)] ) 

2E{lead, sub-lead} 


where 5* is 1 if the tau candidate originates from a jet and 0 otherwise. The tau ID and 
trigger selection criteria for simulated events are modified as follows: the BDT tau ID 
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criteria for the sub-leading tau candidate is removed if the candidate originates from a 
jet, the BDT tau ID criteria for the leading tau candidate and the trigger requirement are 
removed if the leading tau candidate originates from a jet. 

Events in the VE-|-jets control region are selected by a single-muon trigger with a px 
threshold of 36 GeV. The events are required to contain one isolated muon that: has px > 
40 GeV, has < 6 % and is matched to the object that passed the trigger. There 

must be no additional muons or electrons and at least one tau candidate with opposite 
charge to the muon. The remaining contamination from multijet events is suppressed 
by requiring cos A(/)(/r, E™®®) -|-cos A(/)(rhad-vis)< —0.15, which disfavours back-to- 
back topologies where the vector points either in the direction of the muon or the tau 
candidate. The leading-px tau candidate is used to measure the fake-rate. Figures 3(left) 
and 3(right) show 1-prong and 3-prong tau candidates, respectively. The fake- 

rates .RiauXD (including the trigger requirement in the numerator) have a similar behaviour 
but are a factor of two to four lower. The requirement of opposite charge between the 
muon and the tau candidate enhances the contribution of the leading-order qg —)• 
process in which the tau candidate originates from a quark-initiated jet. To evaluate the 
systematic uncertainty from applying these fake-rates to simulated samples with different 
jet origin, the fake-rates are also calculated for events where the tau candidate has the same 
charge sign as the muon. These events have a higher fraction of gluon-initiated jets and 
represent an extreme variation in the jet compostion, resulting in lower fake-rates as shown 
in Fig. 3(left) and Fig. 3(right). A 60% uncertainty is assigned to cover the range of the 
measured fake-rates for events with opposite- or same-sign tau candidates. The uncertainty 
is omitted for lT-|-jets events as they are expected to have the same jet composition as 
events in the control region. The statistical uncertainty from the limited size of the W -|-jets 
control region is also considered. 

7 Systematic uncertainties 

Systematic effects on the contributions of signal and background processes estimated from 
simulation are discussed in this section. These include theoretical uncertainties on the cross 
sections of simulated processes and experimental uncertainties on the trigger, reconstruc¬ 
tion and identification efficiencies; on the energy and momentum scales and resolutions; 
and on the measurement of the integrated luminosity. Uncertainties on the background 
contributions estimated from data are discussed in their respective sections. 

The overall uncertainty on the Z' signal and the Zj^* —)■ eejgplrr background due 
to choice of the PDFs, as, and the renormalisation and factorisation scales is estimated to 
be 14% for a ditau mass of 1750 GeV, dominated by the PDF uncertainty [24]. The uncer¬ 
tainty is evaluated using 90% CL MSTW2008NNLO PDF error sets and also takes into ac¬ 
count potential differences between the following PDFs at the same as: MSTW2008NNLO, 
CTIONNLO, NNPDF2.3 [76], ABMll [77] and HERAPDF1.5 [78]. Additionally, for 
Z/ 7 * —>■ rr, a mass-dependent systematic uncertainty of up to 4% is attributed to elec- 
troweak corrections [24]. This uncertainty is not considered for the signal as it is strongly 
model dependent. An uncertainty of 5% is estimated for diboson production, derived from 
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Figure 3. Tau identification fake-rate measured in VF(—)■ /ri/)-|-jets data events for the BDT loose 
identification working point for (left) 1-prong and (right) 3-prong tau candidates. The fake-rate is 
parameterised in the charge product of the muon and fake tau candidate. Opposite-sign events are 
depicted by black circles and same-sign events by blue stars. The systematic uncertainty covers 
differences due to jet composition and is added to the statistical uncertainty in quadrature. 
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scale, PDF and as variations. A 6% uncertainty on the VF-|-jets normalisation is derived 
from comparisons to data in the lF-|-jets control region used to measure jet-to-tau fake-rates 
in the ThadT'had channel. For tt and single-top-quark production, the uncertainties from vari¬ 
ations in the renormalisation and factorisation scales are in the range of 3-6% [57, 79, 80], 
while those related to the proton PDFs amount to 8% [48, 51, 81-83]. 

The uncertainty on the integrated luminosity is 2.8%. It is derived from a prelimi¬ 
nary calibration of the luminosity scale derived from beam-separation scans performed in 
November 2012, following the same methodology as that detailed in ref. [84]. Comparisons 
of the efficiency of the hadronic tau trigger measured in data and in simulation are used to 
derive an uncertainty of 10% on the trigger efficiency. Differences between data and sim¬ 
ulation in the reconstruction and identification efficiency and the energy scale of hadronic 
tau decays are taken into account. The associated uncertainties for muons and electrons 
are negligible for this analysis. 

The systematic uncertainty on the identification efficiency of hadronic tau decays is 
estimated at low px from data samples enriched in Z —)• rr events, yielding an uncertainty 
of 2-7% depending on the number of tracks and \r]\ of the tau candidate. At high px; there 
are no abundant sources of real hadronic tau decays from which an efficiency measurement 
could be made. Rather, the tau identification is studied in high-px dijet events as a func¬ 
tion of the jet px) which indicates that there is no degradation in the modelling of the 
detector response as a function of the px of tau candidates. Based on the limited precision 
of these studies, an additional uncertainty of 14% • px/TeV for 1-prong tau candidates and 
8% -px/TeV for 3-prong tau candidates is added in quadrature to the low-px uncertainty 
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for candidates with pT > 100 GeV. The reconstruction efficiency for 3-prong tau candi¬ 
dates decreases at high px due to track merging. An uncertainty of 50% -px/TeV above 
Pt = 150 GeV is assigned for 3-prong candidates, derived from data/MG comparisons of 
tracking performance within jets. The energy scale uncertainty for hadronic tau decays 
and jets is evaluated based on the single-hadron response in the calorimeters [66, 69]. In 
addition, the tau energy scale is validated in data samples enriched in Z —)• rr events. 
The systematic uncertainty related to the tau energy scale is a function of r/, px and the 
number of prongs, and is generally near 3%. Energy scale and resolution uncertainties for 
all objects are propagated to the E™'®® calculation. The uncertainty on the due to 

clusters that do not belong to any reconstructed object has a minor effect. 

Table 2 summarises the systematic uncertainties across all channels for the 1750 GeV 
ZggM mass point. In the Thad^had channel the dominant uncertainties on both the signal and 
background come from the tau efficiency and energy scale, while in the riepThad channels 
the statistical uncertainty on the background coming from the fake-factor estimate also 
makes a major contribution. The uncertainties are the same for background and similar 
for the signal for all higher signal mass points, since the same thresholds are used. 
The uncertainties for the lower mass points are typically very similar, except for the tau 
ID efficiency, the 3-prong tau reconstruction efficiency, the Zj^* cross section and the 
statistical uncertainties, which are all a few percent lower, and the uncertainty on the tau 
energy scale for the signal, which can be up to 11% at low mass since the requirement 
is much tighter relative to the Z' mass. The small data-driven uncertainty contribution to 
the signal in the Thad'^had channel comes from jets that are misidentified as hadronic tau 
decays. 


Uncertainty [%] 

Signal 

’^had'Hiad '^fi'^had 

'^e'^had 

Background 

’^had'Hiad '^fi'^had '^e'^had 

Statistical uncertainty 

2.4 

4 

4 

6 

21 

21 

Efficiency 

16 

8 

8 

12 

5 

4 

Energy scale and resolution 

2.9 

5 

5 

10 

11 

9 

Theory cross section 

- 

- 

- 

6 

6 

6 

Luminosity 

2.8 

2.8 

2.8 

2.5 

2.2 

1.9 

Data-driven methods 

0.2 

- 

- 

2.7 

8 

12 

Total 

17 

11 

10 

18 

27 

28 


Table 2. Uncertainties on the estimated .^sgM contribution = 1750 GeV) and the cor¬ 

responding total background contribution in percent for each channel. A dash denotes that the 
uncertainty is not applicable. The statistical uncertainty corresponds to the uncertainty due to the 
limited size of the samples produced via simulation or selected in control regions. The total consists 
of all uncertainties added in quadrature. 


8 Z' signal models 

In this section, the impact on the signal acceptance times efficiency from altering the Z' 
couplings and from including interference between Z' and Zj^* is discussed. The accep- 
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tance times efficiency for a given Z' model is defined as: 

Ae=^^ 

/lint • oB 

where is the expected number of Z' events passing the full analysis selection, uB is 
the Z' production cross section times r+r” branching fraction and C\nt is the integrated 
luminosity. The impact on Ae is presented as a fraction of the SSM value, »4essM- The 
corresponding impact on the acceptance alone. A, is also evaluated by replacing with 
the expected number of Z' events after applying the kinematic selection directly to the 
generated particles before simulation. A Z' model that couples preferentially to third- 
generation fermions is also discussed. 

8.1 Z' signal acceptance 

Changing the fermionic couplings of the Z' from their SSM values can alter the signal 
acceptance of the analysis. Such changes are primarily due to alterations in either the 
tau polarisation or the total Z' decay width. Alteration of the tau polarisation changes 
the tau decay kinematics. Most importantly it affects the visible momentum fraction, 
which enters the analysis through the pT thresholds of the reconstructed visible tau decay 
products and via the threshold on The most extreme impact on the acceptance is 

seen for models that couple only to left-handed or right-handed tau leptons: Z'-^ and 
respectively. Alteration of the quark couplings can impact the acceptance if it alters the 
tau polarisation. However, the maximum impact is much smaller than when altering the 
couplings to tau leptons. As the kinematic limit (due to the collision energy) for high- 
mass Z' production is approached, the signal exhibits an increased fraction of low-mass 
off-shell production. The fraction of off-shell events increases rapidly as a function of the 
decay width. Figure 4 shows Ae for the and Z(^ models, and two models with artificially 
altered decay widths: {T/mz' = 1%) and ZT^g {T/mz' = 20 %), each divided by Ae 

for Zggjyj (T/m^/ ~ 3%). Interference between Z' and Z/ 7 * is not included. The statistical 
uncertainty is typically below 5% but can be up to 14% at low mass. A smoothing is 
applied to reduce fluctuations. For Z^ and the largest impact is observed at low mass, 
where the px and thresholds are much more stringent on the signal. In this case, 
alteration of the tau couplings can lead to changes of up to -|-50% and —25%. The impact 
on the Thad'Thad and riepThad channels are different due to the different effect of polarisation 
on leptonic and hadronic tau decays. For and ZF^^, the impact is most prominent 

at high mass where changes of up to -|- 20 % and —45% are observed. At low mass, Ae only 
changes for widths above 10%. The impact is the same for all channels. For all Z' models, 
the change in A is very similar to that in Ae, indicating that the efficiency is insensitive 
to changes in the Z' couplings. 

The impact of interference between Z' and Z/ 7 * is typically small. For the SSM, 
it leads to a reduction in the expected Z' contribution of up to 10% for mz’ < 2 TeV, 
and up to 35% for the highest mass hypotheses. For Z^, and the impact 

is negligible. For ZF^g the impact can be substantial and is highly dependent on the 
choice of the fermionic couplings. An exhaustive treatment is outside the scope of this 
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article. Reinterpretations of the SSM results for models with large widths should specifically 
calculate the impact from interference. 



Figure 4. Signal acceptance times efficiency for Z^, and divided by the 

acceptance times efficiency for Zggj^ as a function of mz', separately for the ThadThad (solid lines 
with filled markers) and TiepThad (dashed lines with empty markers) channels. The statistical 
uncertainty is typically below 5% but can increase to 14% at low mass. 


8.2 Non-universal G(221) model 

The non-universal G(221) model [8, 20, 21] (also known by other names such as Topflavor) 
is an extension of the SM, containing additional heavy gauge bosons, that 

may couple preferentially to third-generation fermions. The model is motivated by the idea 
that the large mass of the top-quark may suggest that the third fermion generation has a 
dynamical behaviour different from the first two generations. Accordingly, the SM weak 
SU{2) gauge group is split into two parts: one coupling to light fermions (the first two 
generations), 517(2)1 and one coupling to heavy fermions (the third generation), SU{2))^. 
The extended gauge group breaks to the SM 517(2)i+h at a high energy scale, u, and then 
eventually to 17(1)em at the usual electroweak scale, v = 246 GeV: 

SU{2\ X 5[/(2)h X U{1)y 4 5C/(2)i+h x U{1)y 4 C/(1)em • 

The mixing between SU{2)\ and SU{2\ is described by the parameter si4 (/>. The 
and bosons are degenerate in mass; the mass is defined at tree level by sin^ (j) and 

u. Large mixing between r and // leptons has been considered as an additional feature of 
the model, but is ignored here as it would lead to stronger limits via the dielectron and 
dimuon searches. The couples almost exclusively to left-handed fermions, and while 
the coupling strength differs for light and heavy fermions, it is largely insensitive to the 
electric charge or weak isospin, leading to almost universal couplings for all light and heavy 
fermions. 
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Figure 5(left) shows the cross section times t^t~ branching fraction, cj^nU) 
divided by ctBssm- For much of the parameter space crBt^xj is larger than cr;B ssM) peaking at 
moderate values of sin^ cj). For extreme values of sin^ (j) either the cross section is suppressed 
by weakened couplings to light quarks (sin^ i;/) ~ 0) or the branching fraction is suppressed 
by weakened couplings to tau leptons (sin^ 0 ~ 1). Figures 5(middle) and 5(right) show 
the acceptance times efficiency, divided by .4essM) for the Thad'^had and riepThad 

channels, respectively. In general is lower than .AessM- At low mass this is mainly 

due to the left-handed couplings, which result in softer visible tau decays. Near sin^ cj) ^ 0 
and sin^ (/> ~ 1, the acceptance loss comes mainly from the significantly increased decay 
width, which causes a large fraction of the signal to be produced off shell. 



Figure 5. Signal production cross section times r+r“ branching fraction for aB^u, divided 
by crKssM (left) and acceptance times efficiency for AenUj divided by -AessM for the (middle) 
ThadUiad and (right) UepThad channels, as a function of sin^ cji and mz'. 


9 Results and discussion 

A summary of the expected number of events remaining after successively applying each 
selection requirement, up to the threshold, for the signal and dominant background 
processes is given in table 3. Figures 6(left) and 6(right) show the distribution after 
event selection in the Thad'^had and riepThad channels, respectively. The numbers of observed 
and expected events (including their total uncertainties) after applying the thresholds 
in all channels are summarised in table 4. In all cases, the number of observed events is 
consistent with the expected Standard Model background. Therefore, upper limits are set 
on the production of a high-mass resonance decaying to pairs. The acceptance and 

acceptance times efficiency for Zggj^ is shown in figure 7. 

The statistical combination of the channels employs a likelihood function constructed 
as the product of Poisson-distributed random numbers describing the total number of 
events observed in each channel. The probability in each channel is evaluated for the ob¬ 
served number of data events given the signal-plus-background expectation. Systematic 
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Z/ 7 * —>■ TT 

ThadThad channel 

Multijet lV/Z/ 7 *-|-jets 

Top - 1 - diboson 

■^SSM 

Preselection 

276(18) 

611(5) 

64(1) 

24(2) 

10 . 1 ( 2 ) 

OS 

270(18) 

316(4) 

53(1) 

21 ( 2 ) 

9.5(2) 

A(()(ri,r 2 ) > 2.7 

117(2) 

209(3) 

35(1) 

11 ( 2 ) 

9.2(2) 



riepThad channel 




Zj^* —)• TT Jet —> r fake 

Zh* U 

Top - 1 - diboson 

7’ 

^SSM 

Preselection 

46 800(300) 

154 670(130) 

17340(250) 

12 330(70) 

14.3(2) 

OS 

46 300(300) 

111270(120) 

16180(240) 

11830(70) 

13.9(2) 

A(;i(£,T) > 2.7 

32 200(300) 

47650(80) 

12 490(210) 

3530(40) 

13.5(2) 

mx < 50 GeV 

29 490(230) 

22 660(60) 

11240(210) 

808(16) 

8.5(2) 


Table 3. Number of expected signal = 1750 GeV) and background events in the Thad'^had 

and TiepThad channels after successively applying each selection criterion. The statistical uncertainty 
in the least signihcant digit(s) is shown in parentheses. 
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Figure 6. The distribution after event selection in the (left) Thad^had and (right) riepThad 
channels. The estimated contributions from SM processes are stacked and appear in the same order 
as in the legend. The expected contributions from three signals with masses of 750, 1250 

and 1750 GeV are shown, stacked on the total SM expectation. The events observed in data are 
overlaid. The hatched area indicates the uncertainty on the total estimated background. The bins 
have a constant width of (left) 0.153 and (right) 0.184 in log(TOT*). The last bin includes overflow 
events. The inset shows the ratio of the observed events over the total expected SM contribution. 
The statistical uncertainty from the observed events and the expected SM contribution are shown 
on the points and by the yellow band, respectively. The red band depicts the total systematic and 
statistical uncertainties on the SM contribution added in quadrature. 


uncertainties on the expected number of events are incorporated into the likelihood via 
nuisance parameters constrained by Gaussian distributions. Correlations between signal 
and background and across channels are taken into account. A signal-strength parameter 
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niz' 

[GeV] 

'Oiad'Tiad 

Ns Nb N 

TepTiad 

Ns Nb 

N 

500 

400 

1030(170) 

70(8) 56 

570(90) 

49(6) 

42 

625 

450 

650(100) 

40(5) 30 

420(50) 

29(4) 

23 

750 

500 

410(60) 

24.0(30) 18 

270(29) 

18.2(23) 

15 

875 

550 

206(30) 

14.6(20) 11 

152(14) 

11.2(16) 

10 

1000 

600 

119(17) 

9.4(13) 4 

82(8) 

6.7(11) 

6 

1125 

700 

60(9) 

4.0(6) 0 

45(5) 

2.5(4) 

2 

1250 

750 

35(6) 

2.8(5) 0 

27.0(29) 

1.78(32) 

1 

1375 

800 

20.8(34) 

1.93(32) 0 

15.7(16) 

1.24(23) 

1 

1500 

850 

13.4(22) 

1.32(24) 0 

9.6(10) 

0.96(20) 

1 

1625 

850 

8.4(14) 

1.32(24) 0 

6.5(7) 

0.96(20) 

1 

1750 

850 

5.4(9) 

1.32(24) 0 

4.0(4) 

0.96(20) 

1 

1875 

850 

3.6(6) 

1.32(24) 0 

2.70(27) 

0.96(20) 

1 

2000 

850 

2.4(4) 

1.32(24) 0 

1.85(18) 

0.96(20) 

1 

2125 

850 

1.54(28) 

1.32(24) 0 

1.19(12) 

0.96(20) 

1 

2250 

850 

1.02(19) 

1.32(24) 0 

0.81(8) 

0.96(20) 

1 

2375 

850 

0.66(12) 

1.32(24) 0 

0.52(5) 

0.96(20) 

1 

2500 

850 

0.43(8) 

1.32(24) 0 

0.330(34) 

0.96(20) 

1 


Table 4. Number of expected ^ggM signal {Ns), background {Nb) and observed (iV) events in 
the ThadT’had and riepr^ad channels. The signal mass {mz') and corresponding thresholds are 
given in units of GeV. The total uncertainty (statistical and systematic added in quadrature) in 
the least significant digit (s) is shown in parentheses. 


multiplies the expected signal in each channel, for which a positive uniform prior proba¬ 
bility distribution is assumed. Theoretical uncertainties on the signal cross section are not 
included in the calculation of the experimental limit as they are model dependent. 

Bayesian 95% credibility upper limits are set on ctBssm as a function of mz', using 
the Bayesian Analysis Toolkit [85]. Figures 8(left) and 8(right) show the limits for the 
individual channels and for the combination, respectively. The resulting 95% CL lower 
limit on the mass of a decaying to pairs is 2.02 TeV, with an expected limit of 

1.95 TeV. The observed and expected limits in the individual channels are, respectively: 
1.89 and 1.80 TeV (rhad'T'had); 1-59 and 1.59 TeV (r^Thad); and 1.55 and 1.65 TeV (reThad)- 
Alteration of the Z' couplings can impact the signal acceptance as described in section 8.1. 
These changes translate linearly to the limits on uBssm- Limits on the and models 
are shown in figure 8 (right). The impact of the choice of the prior on the signal-strength 
parameter is evaluated by also considering the reference prior [86]. Use of the reference 
prior improves the limit on itBssm by a maximum of 10%, corresponding to an increase of 
20 GeV in the observed mass limit. 

Limits on the non-universal G(221) model are also calculated. The signal contributions 
in the Thad'^had and riepThad channels are rescaled by (tBnu/o''Sssm‘-4£nu/-^£ssm (as derived 
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Figure 7. Acceptance (dashed lines with empty markers) and acceptance times efficiency (solid 
lines with filled markers) for -^sgM a function of the mass. Contributions from the individual 
channels and the full analysis are given. 




Figure 8. Bayesian 95% credibility upper limits on the cross section times ditau branching fraction 
for a Z' in the Sequential Standard Model. The figure shows (left) an overlay of the observed (solid 
lines with filled markers) and expected (dashed lines with empty markers) limits in each channel 
and for the combination, and (right) the combined limit with lu and 2cr uncertainty bands and an 
overlay of the impact of the models. The width of the ^ggM theory band represents the 

theoretical uncertainty from the PDF error set, the choice of PDF as well as the strong coupling 
constant. 


in section 8.2). In addition, the systematic uncertainties are re-evaluated for each point 
in parameter space. Figure 9 shows the region in the parameter space excluded at 
95% credibility: bosons with masses below 1.3-2.1TeV are excluded in the range 
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0.03 < sin^ (j) < 0.5 assuming no /U-r mixing. Indirect limits are also overlaid. 
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Figure 9. Observed 95% CL exclusion in the non-universal G(221) parameter space from the 
combination of the ThadThad and riepThad channels (blue). Indirect limits at 95% CL from fits 
to electroweak precision measurements (EWPT) [27], lepton flavour violation (LFV) [28], CKM 
unitarity [29] and the original Z-pole data [8] are overlaid. 


10 Conclusion 

A search for high-mass ditau resonances was performed using 19.5-20.3 fb“^ of data col¬ 
lected with the ATLAS detector in pp collisions at ^/s = 8 TeV at the LHC. The Thad'^’kad; 
Tfj.Tha.d and TeThad channels are analysed. The observed number of events in regions with 
high total transverse mass are consistent with the SM expectations. Limits are set on the 
cross section times branching fraction for such resonances. The resulting lower limit on 
the mass of a Z' decaying to r+r” in the Sequential Standard Model is 2.02 TeV at 95% 
credibility, in agreement with the expected limit of 1.95 TeV in the absence of a signal. 
The impact on the Z' acceptance from altering the Z' couplings and the total decay width 
is found to be significant, and can be up to 50%. The impact from interference between 
Z' and Z/y* is only important when the total decay width of the Z' is large. Limits on 
the non-universal G(221) model are also established, where Z' bosons with masses below 
1.3-2.1 TeV are excluded in the range 0.03 < sin^ (p < 0.5 assuming no p-r mixing. 
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